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ABSTRACT 

We describe and model high velocity (~ 15,000 km s _1 ) Lya and Ha 
emission from supernova remnant 1987A seen in September and October 1997 
with the Space Telescope Imaging Spectrograph. Part of this emission comes 
from a reverse shock located at ~ 75% of the radius of the inner boundary of 
the inner circumstellar ring and confined within ±30° of the equatorial plane. 
Departure from axisymmetry in the Lya and Ha emission correlates with that 
seen in nonthermal radio emission and reveals an asymmetry in the circumstellar 
gas distribution. We also see diffuse high velocity Lya emission from supernova 
debris inside the reverse shock that may be due to excitation by nonthermal 
particles accelerated by the shock. 

Subject headings: supernova remnants: individual (SNR 1987A) - circumstellar 
matter - shocks 
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1. Introduction 

On May 24, 1997 (day 3743 after core collapse), Sonneborn et al. (1998) observed 
supernova remnant (SNR) 1987A with the Space Telescope Imaging Spectrograph (STIS). 
They discovered very broad (« 20, 000 km s" 1 ) Lya emission in the ultraviolet spectrum. 
Michael et al. (1998, hereafter M98) interpreted this emission as the result of excitation of 
neutral hydrogen atoms in the supernova debris crossing a reverse shock. The observed Lya 
flux was close to the value predicted by Borkowski et al. (1997, hereafter B97) for such a 
model normalized to account for the soft X-ray flux observed by ROSAT. M98 also pointed 
out that the broad Lya emission should be accompanied by Ha emission that should be 
observable both from the ground and with STIS. 

The observations reported by Sonneborn et al. (1998) were taken with a 2" x 2" 
aperture that included the entire inner circumstellar ring of SNR 1987A. With 24 A arcsec 
dispersion, the resulting Lya image seen on the STIS detector is a convolution of the spatial 
(w 1") and velocity (~ 3") structure of the emitting region. M98 found that they could 
fit the observed image with a model in which the shock surface was a prolate ellipsoid 
slightly elongated in the direction normal to the ring plane, but they found evidence for an 
additional source of diffuse emission inside the shock surface. This result was surprising, 
because the model in B97, following Chevalier & Dwarkadas (1995), implies that most of 
the emission should come from a ring shaped surface inside the inner circumstellar ring. 

Here we report a new set of STIS observations of the fast shocks in SNR 1987A, this 
time taken with slits covering only part of the emitting region and including Ha as well 
as Lya. These new observations enable us to distinguish the dynamics and morphology of 
the fast shocks more clearly. We now see that the shock emission is indeed concentrated in 
the equatorial plane, contrary to our interpretation of the earlier STIS observations. We 
also see a diffuse component of emission that may indicate excitation and ionization by 
nonthermal particles. We have not yet detected an unambiguous signal from high velocity 
Nv AA1239, 1243 emission. 

The same set of observations also provides new information on the spectra of the 
supernova ejecta and the circumstellar ring, including a rapidly brightening "hot spot" that 
is apparently caused by the first encounter of the supernova blast wave with the ring. We 
will interpret those data in subsequent papers. 
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2. Observations 

The Supernova Intensive Study (SINS) collaboration obtained STIS spectra of SNR 
1987A in ultraviolet and optical wavelengths in September and October 1997 respectively. 
The ultraviolet spectrum, which covers the wavelength range 1130 - 1720 A, was taken 
on JD 2,450,719.1, or 3869.3 days since explosion, with the G140L low dispersion grating 
(Af « 300 km s _1 at Lya). The resolution along the spatial direction is ~ 0'.'06. Five 
spectral images totaling 11200 sec of exposure were combined. The optical spectrum, which 
covers the wavelength range 5260 - 10260 A, was taken on JD 2,450,727.8, or 3878.0 days 
since explosion, with the G750L low resolution grating (Av ~ 450 km s _1 at Ha). The 
spatial resolution is ~ O'.'l. Two spectral images totaling 4800 sec of exposure were obtained. 
The two images were dithered along the orientation of the slit so that both cosmic ray hits 
and hot pixels on the STIS/CCD were removed simultaneously in the final combined data. 

Figure |l|a illustrates the slit positions on the inner circumstellar ring for the two 
observations. A 0'.'5 slit was used for Lya and a 0'.'2 slit was used for Ha. Both slits were 
offset O'.'l to the SE of the centroid of the ring, with position angles of —147° and —139° for 
the Lya and Ha observations, respectively. 




Fig. 1. — (a) The position of the Lya (wide) and Ha (narrow) slits on the inner ring, N 
denotes the northern side of the remnant, (b) position of the slits on the reverse shock 
surface which lies at ~ 70% of the distance to the inner ring. 

Figure |2|a shows the ultraviolet spectrum in the vicinity of Lya. It displays the 
following components: (1) bright stationary geocoronal Lya emission filling the slit; (2) a 
nearly horizontal streak of intense blue-shifted Lya emission below the centerline but inside 
the lower (NE) side of the ring extending to —15, 000 km s -1 ; (3) a fainter horizontal streak 
of red-shifted emission above the centerline but inside the upper (SW) side of the ring 
extending to 12, 000 km s" 1 ; (4) a vertical band at 4000 - 8000 km s _1 due to diffuse, N v 
AA1239, 1243 emission from nearly stationary N v ions in the vicinity of the ring; and (5) a 
broad component of diffuse emission extending from —15, 000 to 20, 000 km s _1 . 

Figure 0c shows the optical spectrum in the vicinity of Ha. It displays the following 
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components: (1) pairs of bright spots above and below the centerline due to emission by 
the nearly stationary inner ring at [Ol] AA6300,6364, [Nil] AA6548,6584, Ha, and [Sn] 
AA6717, 6731; (2) broad horizontal bands along the centerline due to emission in these 
same lines from atoms in the inner supernova debris excited by radioactivity; and (3) a 
nearly horizontal streak of blue-shifted Ha emission below the centerline but inside the 
circumstellar ring extending to —12,000 km s _1 . 




Velocity [1000 km/s] 




Velocity [1000 km/s] 



Fig. 2.— (a) STIS G140L spectrum of SNR 1987A in the vicinity of Lya. (b) Simulated 
STIS image of Lya emission from the reverse shock surface, (c) STIS G750L spectrum of 
SNR 1987A centered on Ha. (d) Simulated STIS image of line emission from the stationary 
ring and Ha emission from the reverse shock surface. The velocity scales are line of sight 
velocity measured from LMC rest frame. 
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3. Reverse Shock Emission 

According to B97 and M98, we expect to see Lya and Ha emission from Hi atoms in 
the supernova ejecta that are excited by electron and ion impacts as they cross a reverse 
shock. (We do not expect to see Lya and Ha emission from the blast wave because it is 
entering a region where the hydrogen atoms are already nearly fully ionized.) Since the Hi 
atoms are freely streaming, with radial velocity v = r/t, before they cross the shock surface, 
the Doppler shifts of the Lya and Ha are determined by the geometry of the surface. M98 
inferred from the earlier STIS observations that the reverse shock surface was a slightly 
prolate ellipsoid. For such a model, they predicted that a STIS observation of Lya through 
a 0'.'5 slit would appear as an elliptical image (their Fig. Id). 

The actual STIS images of Lya and Ha (Figs. ^|a,c) do not agree well with the prolate 
ellipsoid model for the reverse shock. Instead, they agree better with the model proposed 
by B97, in which the reverse shock emission is confined primarily to an equatorial band. 
We find also that the position and shape of the reverse shock surface vary from the far to 
the near side of the remnant, and probably on smaller scales as well. 

We can fit both the blue-shifted Lya and Ha observations with a model in which the 
shock on the near side has a radius which ranges from 0.73R r i ng on the left side of the (X'5 
slit to 0.70R r i n g on the right side of the slit and extends in latitude from 30° to —20°. 

To fit the red-shifted Lya observations, the shock surface on the far side must have a 
radius 0.76R ring and extend in latitude from 10° to —5°. The red-shifted Ha emission was 
too faint to measure. 

Figure |l]b illustrates the reverse shock surface and its orientation in the slits for the 
Lya and Ha observations. The surface outside of the slits is not constrained by the present 
data and is shown only as an extrapolation. Figures |2|b,d show the resulting model Lya and 
Ha STIS images, respectively. Lines of constant latitude are shown on the simulated Lya 
image assuming that all the emission comes from the center of the slit (see §4). Comparison 
of the model Lya and Ha images demonstrates the advantage of a narrow slit in mapping 
the shock surface. 

When we subtract the simulated images from the actual observations, the streaks 
vanish, leaving residual emission in the Lya observation, which we discuss below in §4. 
The measured flux from the Lya shock emission is (8.6 ± 2.6) x 10" 15 ergs cm" 2 s _1 on 
the blue side and (1.8 ± 0.8) x 1CT 15 ergs cm _2 s~ 1 on the red side. The measured flux 
from the Ha shock emission is (6.0 ± 1.8) x 10~ 16 ergs cm" 2 s" 1 on the blue side and 
< 4 x 10~ 16 ergs cm" 2 s" 1 on the red side. 
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The blue-shifted Lya flux from the shock is greater than the red-shifted flux by a factor 
of ~ 4.8. If we assume that the outer SN debris is spherically symmetric with a power law 
density profile given by p(r,t) oc t~ 3 (r/t)~ 9 (Eastman & Kirshner 1989), we can attribute 
this asymmetry to: the difference in shock radius (which contributes a factor as 1.7); the 
fact that the reverse shock has a greater surface area on the near side (a factor m 2.1), and 
the ~ 0.3 year difference in light-travel time from the near side of the shock surface to the 
far side (a factor 1.3). 

This asymmetry should not come as a great surprise. Gaensler et al. (1997) found that 
the nonthermal radio emission seen inside the inner circumstellar ring was much brighter on 
the northeast side than on the southwest side. The asymmetry in the reverse shock surface 
is not caused by an asymmetry in the supernova debris because then the red-shifted Lya 
emission would be brighter due to the higher mass flux (ram pressure) required to drive 
the shock to a greater radius. The asymmetry therefore represents an asymmetry in the 
circumstellar gas distribution which the ejecta is encountering (the circumstellar gas on the 
far side must be more tenuous than the near side). 

M98 calculated that, for the conditions behind the reverse shock, approximately 1 Lya 
and 0.2 Ha photons will be emitted per Hi atom entering the shock. We can infer this 
emissivity ratio from the observations and our model, which allows us to correct for the fact 
that the Lya observation is taken through a wider slit than the Ha observation. Correcting 
for interstellar extinction factors of 0.15 at Lya and 0.63 at Ha (Fitzpatrick et al. 1990) we 
find, on the blue-shifted side, where the Ha streak is bright enough to measure, this ratio 
to be 4.6 ± 2, consistent with our expectation. 

Since we know the number of Lya photons emitted per H I atom crossing the shock 
surface, we can infer the normalization [A, the density at (r/t) = 10, 000 km s _1 and t — 10 
years] of the SN envelope's density profile from our model for the shock surface. Estimating 
that the reverse shock is propagating outward with velocity V rs = 4, 000 km s -1 (B97) and 
assuming a helium number density n Hc = 0.2nn, we find A = 170 ± 60 amu cm -3 , which 
may be compared with Shigeyama & Nomoto's (1990) model 14E1 (A = 120 amu cm~ 3 ) 
and Woosley's (1988) model 10H (A = 360 amu cm -3 ). If we had taken the value of 
V rs inferred from radio observations of the remnant (Gaensler et al. 1997), A would be 
decreased by « 10%. 
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4. Diffuse Emission 



When we subtract the model images from the actual images, a diffuse component of 
high velocity Lya emission remains. This component was evident in the earlier slitless 
observation by Sonneborn et al. (1998), and led M98 to interpret those data with a model 
in which the reverse shock surface had the shape of a slightly prolate ellipsoid. Now, with 
the present observations, we can see clearly that the Lya and Ha emission from the reverse 
shock is confined to the equatorial plane, and that the residual Lya emission comes from 
a volume interior to that surface. A corresponding diffuse component of Ha presumably 
exists, but it is too faint to be seen with the present exposure. 

To analyze the distribution of the diffuse Lya emission, we assume that the emitting 
gas is everywhere in free expansion. If so, and if the slit were very narrow, there would be a 
unique mapping from the STIS image to depth within the supernova debris. From such a 
mapping, we could reconstruct a two-dimensional image of the supernova debris intercepted 
by the slit. Figures |3|a-c are just such images, where the coordinates are radius and latitude 
measured from the equatorial plane of the inner circumstellar ring. The curves of constant 
latitude are not radial lines because the slit is off-center. 




Fig. 3. — Remapping of Lya emission where appropriate, see text (§4). (a) observed 
emission, (b) simulated reverse shock emission, and (c) interior blue-shifted emission (panel 
(b) subtracted from panel (a)). 



Figure |^a represents the Lya data of Figure 0a in such a coordinate system, in which 
we assume that the slit is a line centered in the actual slit. Since the actual slit has width 
0'.'5, the mapping is not unique; this introduces a "point-spread function" for the mapping 
represented approximately by the error circles in Figures §a-c. 



Since we cannot distinguish the contribution of red-shifted Lya from that of Nv 
A A 1239, 1243 with the present data, we cannot interpret the signal on the lower right hand 
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side of Figure [3]a. 

Figure |3|b represents the model for Lya emission from the reverse shock surface 
corresponding to Figure [l]b, including the broadening due to the 0'.'5 slit. In this coordinate 
system, we can see clearly that the blue-shifted Lya emission comes from a reverse shock 
surface that lies primarily above the equatorial plane. 

Figure 0c shows the residual flux on the blue side of Lya after subtracting the 
contribution from the reverse shock from the data. This diffuse Lya emission comes from 
supernova debris inside the reverse shock surface and appears to have a roughly uniform 
radial distribution. Like the Lya emission from the reverse shock surface, the diffuse Lya 
emission is strongest within latitudes 0° to 40°; but we see significant diffuse emission from 
-20° to 90°. 

The flux of diffuse blue-shifted Lya seen through the Of! 5 slit is 
(5.7 ± .3) x 10 -14 ergs cm _2 s _1 , or w 7 times the flux of blue-shifted Lya from the 
reverse shock. Each pixel in the spectrum observes a volume element in the debris so 
we are able to measure the volume emissivity of the diffuse interior emission. Correcting 
for interstellar extinction and assuming a distance of 50 kpc to the remnant we find this 
emissivity to be « 4 x 10 -19 ergs cm -3 s _1 . 

What causes this diffuse emission? Since it is concentrated in the same latitude range 
as the reverse shock emission, one is tempted to look at the reverse shock itself as the 
source of the excitation. We can rule out one obvious candidate — soft X-rays emitted by 
the shock — because the soft X-ray luminosity observed by ROSAT (Hasinger et al. 1996) 
is roughly an order of magnitude too faint to account for the diffuse Lya. Likewise, the hot 
spot is still too faint to cause significant diffuse Lyct emission. We suspect that the diffuse 
Lya emission is caused by nonthermal particles that are accelerated at the reverse shock 
and diffuse into the freely expanding supernova envelope. We know that such particles 
are present from observations of nonthermal radio emission (Gaensler et al. 1997). The 
reverse shock may be a fairly efficient source of nonthermal particles (Jones & Ellison 1991). 
Moreover, those nonthermal particles that diffuse into the supernova envelope and deposit 
their energy may produce Lya photons with efficiency approaching « 30%. 

5. Future Investigations 

Further observations with STIS will give us a much clearer picture of the complex 
dynamics and kinetics of this extraordinary event. It would be valuable to observe the 
system again in the same slit configurations to measure the rate of brightening in Lyct and 
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Ha emission from the shock surface. We should map the high velocity emission from the 
entire system, especially near the east side where the radio lobes are brightest. 

As M98 have emphasized, observations of the line profile of high velocity Nv 
AA1239, 1243 emission will provide a unique probe of the kinetics of a collisionless shock. 
We believe that the emission seen on the red side of Lya contains a significant contribution 
from N v AA1239, 1243, but we cannot distinguish this emission from Lya emission in the 
present observations. To do that, we must observe the high velocity emission from the same 
part of SNR 1987A with STIS in both Lya and Ha. The exposure should be deep enough 
to see the diffuse emission clearly in Ha as well as Lya, and the slit should be narrow 
in order to separate the spatial structure from the velocity structure and to isolate the 
N v AA1239, 1243 emission from stationary gas. Then, since the ratio of Lya/Ha should 
be constant, we would be able to subtract the Ha image from the Lya image, leaving a 
remainder that is due to N v AA1239, 1243. 
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